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ABSTRACT 

There are two purposes of this project. One is to 
determine whether linear programming techniques can 
improve the performance in handling design optimiza- 
tion problems with a large number of design variables 
and constraints relative to the feasible directions algo- 
rithm. The second purpose is to determine whether 
using the Kreisselmeier-Steinhauser(KS) 1 function to 
replace the constraints with one constraint will reduce 
the cost of the total optimization. Using the software 
program, CONMIN 2 , reference cases are run with both 
the linear and non-linear options. Next the same test 
case is run using the linear programming subroutine, 
LINPR1 3 * from the math library. Comparisons are then 
made between the solutions obtained from both sub- 
routines, CONMIN and LINPR1. 

PROCEDURE 

A simple problem of a hub with 12 spokes was used 
as the test example (see Appendix). This problem has 
12 design variables and 24 constraints. The calculations 
were done on a DEC MICRO VAX II workstation 
using code written in FORTRAN 77. 

Using CONMIN, results were obtained using the 
non-linear and linear options. Next the 24 constraints 
were replaced by one, a KS-based cumulative con- 
straint, and results obtained, again for both the non- 
linear and linear options (figure 1). 

Since the value of rho in the KS function influences 
the result, different values of rho were used in the KS 
function. The comparisons between the linear and non- 
linear CONMIN solutions using the 24 constraints and 
one constraint(when using the KS function) are in 
figure 2. 

After obtaining the results using CONMIN, the 
non-liner problem was turned into a linear program- 
ming one to be solved using the linear programming 
subroutine LINPR1. This is a library routine from the 
math library and uses the simplex method. The results 
using this routine are compared to the CONMIN 
results in figure 3. 

CONCLUSIONS 

The optimal value of the objective function was 
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Figure 1 CONMIN, linear, non-linear for 
normal and KS 
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Figure 2 CONMIN, non-linear, linear, KS 

practically the same when using CONMIN with or 
without linearization and with or without the KS func- 
tion as the cumulative constraint (figure 1). When the 
KS cumulative constraint was used, the optimal objec- 
tive function was influenced by the rho factor (figure 2). 
The CONMIN optimized objective function was 
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consistently lower than the one obtained from the linear 
programming routine, regardless of the use of the KS 
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Figure 3 CONMIN vs. LIN PR 1 

function (figure 3). This is an unexpected and impor- 
tant finding of this study. 

In terms of efficiency, all the runs were comparable 
in the number of function evaluations needed. 
However, there is a reduction of memory required by 
CONMIN when the constraints are replaced by a single 
cumulative constraint using the KS function. 
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APPENDIX 

Problem Formulation: 

y 


X 


Nomenclature 


NS 

number of spokes 


NLC 

number of loading cases 


Oi+l 

= Oi + 2rc/NS 


E 

Youngs modulus 


Ai 

cross-sectional area of rod i 


R 

radius of the circle = length of each rod 

uJ* 

displacement of the hub along x 
loading case 

for jth 

“V 

P'lJP'y 

displacement of the hub along y 
loading case 

for jth 

load components along x, y for jth 

loading 


case 


bi 

=EA i/R 


CTai 

allowable tension stress 


CTac 

allowable compression stress 


Details 

of the Analysis 



NS 



*11 — y ^bjCQ^ai 
% 

(1) 


NS l 



*12 = J2 2 binn2at 

ft 

(2) 


NS 



* 2 2 = y bjsin2oti 

(3) 


Displacements for loading casesj 



DET = knk22-k 2 i2 
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(4) 




where 


iPx = (P’fc - P , yki2)/DET 

(5) 

u^ y = (P’yku - P*Jci 2 )/DET 

(6) 

Strain in rod i for loading case j 


e^i = (-u^cosai - i^ysinc^i/R 
Stress in rod i for loading case j 

(7) 

tu 

ii 

~is) 

(8) 

There are NS* NLC stresses c^i 
Material volume 



9) 


DETAILS OF THE OPTIMIZATION USING CONMIN 

min V(Ai) (10) 

Ai 

subject to the constraints 

gm= (cHi/Oa - 1) <0 (11) 

m = 1 — > (NS» NLC) (12) 

where ct» = a»t. of cr 1 ; > 0, otherwise cr a = cr ac 
Amin 5 Ai 

Numerical Data for Test Cast 

NS = 12 

NLC = 2 

Ai = 1cm 2 

R = 100cm 

E = 20* 10 6 N/cm 2 

loading case 1 

p‘x = 40000 N; P ! y = 0 

loading case 2 

p 2 x= 40000N; Pjr = o 

Oat = 80000 n/cm , cr ac = -40000 N/cm 2 

Amin = .1cm 2 

Using the K-S function In the Optimization. 

The optimization is repeated using the cumulative K-S 
function in place of equation 11. 
gm = KS 


KS = g m *x + -In 

P 


■ NSNLC 
. i=l 


(13) 


DETAILS OF THE OPTIMIZATION USING LINPR1 
Turning the Optimization Problem into a Linear 


Programming One 

Introduce a new variable 

Xi = 1/Ai (14) 

Compute derivatives 

dV/dXi and 3gWdXi (15) 

Because of 14, equation 15 becomes 

av/axi= av/3Ai (-i/x 2 o ^ av/aA, (-a 2 o (i6) 

and dgrJdXi = agim/aA.t-A 2 ,) (17) 

Using equation 9 put 3V/3Aj = R in equation 16 

then 

av/axi= r(-a 2 o (is) 


Approximate Linear Optimization Problem 

Let V°, g° m be the values at the initial Xi = X°i 
Approximate V(Xi), gm(Xi) by extrapolation 


v = v°+ av/axi(Xi - x°o (i9) 

gm = g°m + agm/aX,(Xi - X°i) (20) 

The approximate problem is: 

min(V° + av/aXi(Xi-X°0) - (21) 

Xi 

STOC 

g° m + agm/ax,(Xi - x°o < o 

where (22) 


pXo < (Xi - X 10 ) < (1 + P)X°i (23) 


3 


equation 23 is a move limit that does not allow Xi 
to move too far from X°i . Initially p = .2. 
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